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We present numerical studies of first-order and continuous filling transitions, in wedges of arbitrary 
opening angle tp, using a microscopic fundamental measure density functional model with short- 
ranged fluid-fluid forces and long-ranged wall-fluid forces. In this system the wetting transition 
characteristic of the planar wall-fluid interface is always first-order regardless of the strength of the 
wall-fluid potential In the wedge geometry however the order of the filling transition depends not 

only on £„ but also the opening angle xp. In particular we show that even if the wetting transition is 
strongly first-order the filling transition is continuous for sufficient acute wedges. We show further 
that the change in the order of the transition occurs via a tricritical point as opposed to a critical- 
end point. These results extend previous effective Hamiltonian predictions which were limited only 
to shallow wedges. 


I. INTRODUCTION 

Fluid interfacial phenomena have been a topic of intense research within statistical physics for over three decades - 
see for example the review articleJ^®. It is now well established that confining a fluid or, equivalently, placing it in a 
strong external potential may induce novel phase transitions, scaling behaviour and criticalit y distinct from those of the 
bulk fluid. Well studied ex ampl es of this include wetting transitions at planar wall d^^ l ^b oUlJ and capillary condensation 
in parallel plate geometrie d^^bs i and capillary grooveaSHlZI Por wetting phenomena it was quickly established that the 
order of the transition depends very sensitively on the range and strength of the competing wall-fluid and fluid-fluid 
forces. In particular, in three dimensions, continuous (critical) wetting transitions require, in general, a fine tuning of 
the interaction strengths, meaning that in experiments and in model calculations alike, the transition is most often 
first-order in character. On non-planar substrates, however, geometrical effects are also important, and may strongly 
influence the order of any phase transition. This is well illustrated by the filling transition occurring for a fluid confined 
in a linear wedge formed by two planar walls meeting at an opening angle ip- In many respects the filling transition 
is a missing link between wetting and capillary condensation, connecting these apparently two distinct phenomena, 
as well as showing several novel features. Macroscopic arguments dictate, that a wedge in contact with a bulk vapour 
at two phase coexistence, is completely filled by liquid if the contact angle 0 <0/ satisfieJi^II^ 

0f(T) = ^. ( 1 ) 

where (tt — ip)/2 is often referred to as the wedge tilt angle a. The filling transition corresponds to the divergence in 
the adsorption as 9{T) — a —)■ O’*" and can be induced either at fixed T by increasing the tilt angle (i.e. making the 
wedge more acute), or by increasing T ^ Tf causing the contact angle to decrease until the condition 0{Tf) = a is 
fulfilled. Wedge filling therefore precedes any wetting transition in the sense that Tf <Tyj where is the wetting 
temperature at which 6(Tt„) = 0. Indeed wedge filling does not actually require there to be any underlying substrate 
wetting transition since it needs only the familiar phenomena of partial wetting. This makes filling transitions easier 
to observe than wetting since one need only tune the substrate geometry rather than the details of the intermolecular 
interactions required to make the contact angle vanish. On the other hand, assessing the order of the filling transition 
is more difficult than for wetting since there is no analogue of the macroscopic contact angle whose measurement 
would indicate the order of the phase transition. However, the order of filling transitions can be distinct from that 
of wetting and is key to understanding the more subtle aspects of the phase transition. These include strongly 
enhanced interfacial fluctuations and, in two dimensions, hidden connections with critical wetting referred to as 
wedge covarianc^^lHSS] these predictions arose initially from studies based on very simple effective Hamiltonian 

models which generalized the standard Capillary-Wave analysis of wetting transitions to the wedge geometry. More 
recently however these predictions have been tested using microscopic models both at mean-field level and beyond. 
For example in two di mens ions, in addition to the known solution for wedge filling in the square lattice Ising model 
at a right-angle corneJ^^Hl^ the transition has been studied within a field theoretical continuum model of fluid phase 
coexistence, which admits an exact solution for arbitrary opening angle^^. This has shown that the wedge covariance 
has a deeper relation to the Lorentz invariance of quantum field models in 1 -|- 1 dimensions. In addition in three 
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dimensions it has been possible to test predictions for the order of filling transitions in simulationj ^^^^ and using 
square-gradient theorjl^ and modern density functional models based on fundamental measure theorji2212U This has 
illustrated, that for right angle wedges the filling transition may indeed be continuous even though the underlying 
wetting transition is first-order. Interestingly the mechanism for this change in order appears to be even more general 
than originally thought based on simple effective interfacial models. 
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FIG. 1: Schematic illustrations of geometrical cross-sections for a) planar wall, b) rectangular wedge, c) open wedge, d) acute 
wedge. Translation invariance is assumed along the y-axis. 

The purpose of the present work is to extend our recent density functional studies of wedge filling at right angle 
corners to more general opening angles. In particular we wish to show that filling transitions that were observed to 
be first-order (continuous) for a right angle wedge, may become continuous (first-order) by making the opening angle 
smaller (larger). Our article is arranged as follows: In the next section we present our microscopic model and the 
derivation of a long-ranged external potential arising from dispersion-like forces, for a three dimensional wedge with 
arbitrary opening angle. We then consider a right-angle corner with two different strengths of the wall-fluid potential, 
which give first-order and continuous filling respectively. These are then studied for different opening angles showing 
the change in order as the wedge is made more acute and open, respectively. We conclude our paper with a discussion 
of the mechanism regarding the change in the order of the transition. 


II. DENSITY FUNCTIONAL THEORY 

Within classical density functional theory the equilibrium density profile p(r) is obtained from minimization of the 
grand potential functional 

^[P] =-^[p]+y dr(U(r)-/r)p(r). (2) 

where V (r) is the external field and /r is the chemical potential. All the information about the fluid model is contained 
in the intrinsic free energy functional IF[p{y)\ which is often split into an ideal gas and excess contribution. Thus, 

^[P] = “^id[p] + -^ex[p] , (3) 

where Jid[p] = ksT f drp(r) [ln(A^p(r)) — l] and A is the thermal de Broglie wavelength that can be set to unity 
without loss of generality. 

In the spirit of van der Waals theory, the excess term is treated in a perturbative manner, and is separated into a) 
a contribution modelling the repulsive hard-sphere (hs) core and b) a contribution from the attractive part u(r) of 
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the fluid-fluid intermolecular potential which is treated in simple mean-field fashion. Hence we write 

J'exM = ^ J drp(r) J drV(r')u(|r-r'|). 

where, in our analysis, u(r) is taken to be a truncated Lennard-Jones-like potential 

( 0; r < cr, 

Mr) = < -4e(y)® ; a <r <ra, 

[ 0; r > r^. 


(4) 


(5) 


which is cut-off at Tc = 2.5 cr. The hard-sphere term J-hs[/o] describes the repulsion between the fluid particles of 
diameter cr which is approximated using Rosenfeld’s fundamental measure theory aJ^ 



dr$({na}) ■ 


( 6 ) 


where the {n„} are six weighted densities. 

In general, the external potential V (r) for arbitrary wall shapes can be constructed by integrating a two-body 
wall-fluid potential 4’wir) over the volume V of the wall which is assumed to be of uniform density 


H(r)=pu, [ dr'(/)u,(|r - r'l), 

Jv 


where in our study (j)w{r) is taken to be 


(l^wir) = -4eu, 



r > a. 


(7) 

( 8 ) 


In addition we impose a hard wall repulsion V{r) = oo whenever the distance from surface of the wall is less than 
cr. We work always in three dimensions but assume translational invariance along the y axis, so that the potential is 
only a function of the Cartesian coordinates x and z (see illustrations in Fig. I). 

For the simplest case of a planar wall occupying a half space z < 0 (see Fig. la), the potential reduces to a pure 
one-dimensional power-law 


V„{x,z) 



z > a 


(9) 


where = —^'xSwPw'X^ measures the strength of the interaction. Similarly, for a right angle corner (Fig. lb) the 
potential maybe writterP^l 


V^l2{x,z) = ayj 



2z* + + 2x‘^ 1 


x,z > (7 


( 10 ) 


where the subscript refers to the opening angle '0. Notice that far from the apex, x —>■ oo or z —>■ oo, this potential 
reduces to a pure power-law characteristic of the planar wall (|^. The potential (10) was used in our previous studies 
of wedge flllin^2213il Here, we extend this analysis to more general wedges whose potential t0/,(x,z) can be readily 
obtained by integrating the pair potential 4>w(r) over a volume of triangular cross-section which is added either to 
14(x, z) or to 14 / 2 ( 2 :, z). For acute wedges, with 0 < 7r/2 (see Fig. Ic), the attractive part of the potential is 


t4(x, z) = Uu 


1 


cosec^'f/: 

(x — Z cot '0)3 


6x^z^ cot^ 0 -I- 3z^ cot^ 0 — 6x^z cot 0 -I- 2z^ + x^z^ + 2x 


41 


2(x — z cot'ip)^z^\M~+M 


( 11 ) 


with the hard wall repulsion applying within the wall and a distance cr from it. It is clear that V^{x,z) reduces 
immediately to the expression of 14 / 2 ( 2 :, 2 ) for 0 = 7r/2. 

For open wedges (see Fig. Id), corresponding to 0 > ■7r/2, on the other hand, the potential is more conveniently 
written as 


14 (x, z) — Oiiu 


6x^z^ tan^ 0 -I- Sx'^ tan^ 0 — 6z^x tan 0 -|- 2x"‘ -I- x^z^ -|- 2z"‘ 
2(z — X tan 0)3x3-0x24-4^ 


sec3 0 2x^ -I- x^z^ -I- 2z^ 1 

(z - xtan0)3 2x^z^\/xMM 

( 12 ) 
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together with the appropriate hard wall restriction. It is easy to verify that this recovers the potential for the right 
angle corner when '0 = 7r/2 and also the planar wall when ip = n. 

The grand potential functional n[p] is minimized numerically on a two dimensional Cartesian square mesh of grid 
size 0 . 1(7 with appropriate boundary conditions. We first determine the equilibrium profile for a planar wall Ptt{z) 
at temperature T and chemical potential /i. This one dimensional density profile is than imposed as a boundary 
condition on the two dimensional density along the normals ni and 112 at a distance L = 40 cr from the apex along 
each wall. Previous studies of filling at a right angle corner have shown that distance L is large enough to avoid 
significant finite size effects and mimic the interface between the wall and the reservoir fluid. 



FIG. 2: Illustration of the finite-size domain used in the numerical minimization of the grand potential functional in acute and 
open wedges. Along the normals ni and 112 the density is fixed to that of the planar wall-fluid interface, P-k{z), to mimic the 
interface with a bulk vapour. The distance L is set at L = 40a which is much larger than the wetting film thickness at a planar 
wall. 


III. RESULTS 

We work at two phase bulk coexistence and at subcritical temperatures T < Tc which, for our truncated Lennard- 
Jones-like potential, occurs at ksTc/e = 1.41. Throughout this work we consider two wall-fluid interaction strengths 
corresponding to i) = 0.8e and ii) e^, = 0.9e. We begin by considering the wetting properties of each planar 
wall-fluid interface determining the temperature dependence of the contact angle d[T), the wetting temperature 
and, from the numerically determined binding potential, the order of the wetting transition. For both values of the 
transition is unequivocally first-order. This is to be expected since the wall-fluid interaction is long-ranged while the 
fluid-fluid potential is effectively short-ranged. We then consider a right-angle wedge and, from determination of the 
free-energy and adsorption, locate the filling transition temperatures Tf which are shown to be completely consistent 
with the thermodynamic result Q. However the filling transitions are now of different order; the transition for the 
weaker potential, for which and Tf are closer to is continuous, in contrast to the stronger potential for which 
the transition is first-order. We then investigate what happens to the location and order of these transitions when the 
opening angle ip is varied. In particular for each value of we determine the value of the opening angle at which the 
order of the filling transition changes. Finally, from a numerically constructed wedge binding potential we are able to 
determine whether the change in order occurs via a tricritical or critical-end point. 


A. The planar wall-fluid interface 

We consider first a planar wall of infinite area A, with the power-law potential (|^, in contact with a bulk fluid of 
volume V at pressure p. Imposing that p(oo) = pi or p(oo) = pg (corresponding to the bulk liquid and gas densities 
respectively) determines the equilibrium profiles p{z) and surface tensions 7 = (H + pV)/A for the wall-liquid (wg) 
and wall-gas (wg) interfaces respectively. Then from Young’s equation = 7 mi+cos 07 ;g, where 7 /g is the liquid-gas 
interfacial tension, we determine the temperature dependence of the macroscopic contact angle 0(T). This is shown 
in Fig. and indicates the presence of a wetting transition occurring at fct,Tu, = 1.37e for e^, = 0.9e and = lAe 
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FIG. 3: Temperature dependence of the contact angle for the wall strengths = O.Qe and = 0.8e showing wetting transitions 
at fcfjTu, = 1.37eand kbT^ = lAe for respectively. In the inset is shown a log-log plot illustrating the vanishing of 9 in the 
vicinity of each T^. This is consistent with the expected first-order singularity 9{T) ~ {T^ — T)® - the straight lines have slope 
equal to 1/2. From the temperature dependence of the contact angle one may also read off the macroscopic prediction for the 
location of the filling temperature in wedges of opening angle ip according to 9(Tf ) = (tt — ip) 12. 


for Ew = O.Se. The asymptotic behaviour as T — )■ is consistent with 9{T) ~ (T^, — T)^, indicating that the wetting 
transitions are both hrst-order, as expected (see inset). As a further check on this we have numerically determined 
the binding potential W{i) = — "fu]i — "fig, where n{£) is the grand potential of a wetting layer constrained to 

be of thickness i (see Fig. |^. These show an activation barrier between the partially wet and completely wet states 
close to Tu, confirming the first-order nature of the transitions. More generally, from the plots for 0{T) we can now 
test the thermodynamic prediction 0 for the location of the filling transition in wedges with different opening angles, 
beginning with the right angle corner. 


B. The right angle wedge 

We now consider the filling transitions in a right-angle wedge, corresponding to the potential z), for the two 

different values of ew Using the plots of 0{T), the macroscopic result Q, predicts that fesTy/e « 1.33 for Su, = O.Qe 
and ksTf/e « 1.375 for = 0.8e. First consider the excess adsorption, defined for the wedge-gas interface, by 

r = ^ J dx J dz [p{x,z) - pg] . (13) 

For large enough coverage, i.e., for mesoscopicaly large values of the excess adsorption, F oc £^(pi — Pg)/L'^ where 
is the perpendicular distance of the liquid-vapour interface from the wedge apex, the stronger wall-fluid potential, 
Ew = 0.9 the adsorption shows a jump, close to the predicted value oi Tf, between two states corresponding to 
small and large coverage of liquid near the wedge apex. This is reflected in the temperature dependence of the grand- 
potential which shows that two branches cross at Tf = 1.335 which is in excellent agreement with Q (see Fig.j^. The 
coexisting prohles are shown in Fig.|^ Note that in the higher coverage state the interface is essentially flat and meets 
the walls at the correct contact angle 9{Tf) « 7r/4. The coverage of this state scales with corresponding, in the 
thermodynamic limit, to a completely filled wedge. For the weaker wall potential however there is no coexistence and 
the adsorption (and grand potential) has a single branch which rapidly, but smoothly, increases as the temperature is 
raised towards Tf (see Fig. [^. This indicates that the filling transition is continuous despite the fact that the wetting 
transition is first-order. 
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FIG. 4: Numerically determined binding potential functions for the planar wall-fluid interface. The upper panel (a) refers to 
wall strength Sw = O.Qe and T < T^j, while the bottom panel (b) refers to e-w = O.Se and T > Tw The presence of the potential 
barrier indicates that the wetting transition is of hrst-order in both cases but note that the barrier is an order of magnitude 
smaller for = 0.8e consistent with T™ being much closer to Tc- 


C. Non-rectangular wedges 


Consider first the stronger wall, = 0.9^, for which the filling transition is hrst-order for the right-angle wedge, 
tp = 7r/2. We now close the wedge so that the opening angle = 45 and repeat our analysis, determining the 
temperature dependence of the adsorption and grand potential after minimizing S4[p] from low-density and high- 
density initial conhgurations. Once again we hnd hysteresis in the adsorption but this is now much diminished 
indicating the transition is only weakly hrst-order. Correspondingly there are two branches to the grand potential but 
these now meet almost tangentially at ksTj/e = 1.28. Again this value of Tf is in full agreement with the macroscopic 
prediction 0. The coexisting states at the hlling temperature are displayed in Fig. Note that the small decrease 
in the adsorption with increasing T in the high coverage phase is simply due to decrease in the pi — pg and does not 
rehect the position iyj of the interface which is saturated. Decreasing the opening angle further reduces the hysteresis 
which eventually vanishes when the wedge is very acute, when ip « 20. This is illustrated in Fig. which shows only 
a smooth but rapid increase of the adsorption as T —>■ where ksTf/e « 1.2. In this highly acute wedge the filling 
transition has therefore become continuous. 

The opposite happens for the weaker potential = 0.8 for which the filling transition was already continuous for 
the right angle wedge. In this case, opening the wedge eventually induces hysteresis when ip « 120° (see Fig. 11) 
indicating that the filling transition, like the underlying wetting transition, is first-order. For ip = 120° the crossing of 
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FIG. 5: Adsorption (a) and excess grand potential density (b) as a function of temperature for a rectangular wedge with 
wall strength e™ = 0.9e. Two distinct branches are obtained by minimization of from initial low-density high-density 
configurations respectively. The branches in the grand potential cross at ksT = 1.335 e consistent with the thermodynamic 
prediction obtained from Eqn. (1). The hysteresis is indicative of a first-order filling transition. 


the branches of the grand potential determines ksTj = 1.385£ which is in perfect agreement with Q. The coexisting 
states at Tf are shown in Fig. 12 


IV. DISCUSSION AND CONCLUDING REMARKS 

In this work, we have presented a non-local density functional study of filling transitions in open and acute wedges, 
extending previous studies which were restricted to right-angle corners. In our model the walls of the wedge themselves 
exhibit a wetting transition (at temperature T^) which is always first-order in nature regardless of the strength of the 
wall-fluid interaction . We have found that in the wedge geometry, the location of the filling transition temperature 
Tf is always in agreement with the thermodynamic prediction, 0{Tf) = {'K — ip)/2, indicating that Tf can be arbitrarily 
lowered below T^, by decreasing the opening angle ij}. In addition we found that by reducing the opening angle one 
can always drive the filling transition second-order implying that the adsorption continuously changes from micro- to 
macroscopic ai Tf. This generalizes our earlier studies of filling at right angle corners and shows that the change in 
order is not restricted to transitions in the proximity of Tc- 

Our central finding, that it is possible to induce continuous filling for sufficiently acute wedges, is in partial agreement 
with longstanding predictions of simple effective Hamiltonian theory for filling in open wedges (with -0 « TrJPSl. Such 
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FIG. 6: Coexisting equilibrium density profiles for a rectangular wedge with the wall strength = 0.9e at the filling temperature 
ksTf = 1.335e. The low-density state (a) corresponds to the adsorption (solid) lines in Fig. 5 and the high density state (b) 
corresponds to the desorption (dashed) lines in Fig. 5. In the high-density state, the liquid vapour interface meets the wall at 
the contact angle 6 = 7r/4 in agreement with the macroscopic prediction of Eq. 0 - 
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FIG. 7: Temperature dependence of the adsorption in a rectangular wedge with wall strength = 0.8e. The curve has its 
maximum gradient at ksT/e ~ 1.38 very close to the location of the filling transition ksTf ~ 1.378e predicted by Eqn. 0 - 
There is no hysteresis indicating the filling transition is continuous. 


interfacial models also predict a change in order from first to continuous filling, when the filling temperature Tf is 
sufficiently below that there is no activation barrier in the binding potential W{€}, defined for wetting at the planar 
wall. Within such models this means that for wedges made from walls that only exhibit first-order wetting, the filling 
transition is first-order when the wedge is open, but continuous when it is sufficiently acute. This is in qualitative 
agreement with our findings. However the mechanism behind the change of order cannot be exactly the same as that 
within the interfacial Hamiltonian description. This is because within our present study only the wall-fluid potential 
is long-ranged which means that binding potential for wetting at a planar wall always has an activation barrier. The 
activation barrier would not be present if one could induce a change in the sign of the leading order term in the 
binding potential (Hamaker constant) which in turn needs a balance between the strengths of long-ranged wall-fluid 
and long-ranged fluid-fluid potentials. Therefore, if one strictly applied the predictions of the interfacial Hamiltonian 
model to the present system then the filling transition would always be first-order albeit very weak since the strength 
of the activation barrier is rather small at Tf. 

So what is the reason for this discrepancy? One option is that the finite-size restrictions in the present numerical 
study have rounded the filling phase transition, which for larger domain sizes L would be (weakly) first-order. While 
we cannot completely rule out this possibility, it is notable that the numerically determined location of the filling 
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FIG. 8: Adsorption (a) and excess grand potential density (b) as a function of temperature for a wedge with opening angle 
^|) = 45° and wall strength = 0.9e. The filling transition is located at ksTf « 1.28£ and is weakly first-order. 

transition is always in excellent agreement with the thermodynamic result 0{Tf) = a. The location of the transition 
therefore is certainly not strongly influenced by the finite-size. As a check on this we have repeated our analysis of 
filling a rectangular corner for Sw = 0.8 e for the much large domain size up to L = lOOcr which again shows only a 
smooth increase in the adsorption consistent with continuous filling. The second option, which appears more likely 
to us, is that the original effective Hamiltonian description does not capture all the details of the filling transition. 
There are indeed plausible reasons for believing this since the original interfacial Hamiltonian model is only applicable 
to shallow wedges and to filling temperatures Tf far from Tc where a simple sharp-kink description of the interface 
structure is reliable. If the wedge is very acute or if Tf « then a sharp-kink approximation ceases to be valid. 
However these are precisely the conditions where we find a change in the order of the filling transition. Extending 
the original effective Hamiltonian model of filling to these regimes requires, at least, both a soft-kink treatment of 
the non-planar interface and a fully non-local description of the interface-wall potential. In addition in acute wedges, 
packing and volume exclusion effects are almost certain to play an important role and are clearly visible in the density 
profile. For continuous and weakly first-order filling transitions, the free-energy landscape, determining for example 
the energy cost of maintaining a coverage of order F oc is so shallow, that any extra stress on the liquid-vapour 
interface may strongly effect the phase transition. Given that packing effects are completely neglected in the original 
effective Hamiltonian theory it appears to us highly likely that this is the source of any new physics within the 
microscopic density functional description of filling. Incorporating all these features into the interfacial Hamiltonian 
theory is extremely challenging, indeed so much so, that a microscopic density functional treatment is a much more 
tractable way of studying the problem. Our results suggest that further work is required to understand how packing 
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FIG. 9: (Color online) Coexisting equilibrium density profiles for a wedge with opening angle ip = 45° and wall strength 
Ew = 0.9£ at the filling temperature ksTf = 1.28e. 
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FIG. 10: Temperature dependence of the adsorption for a wedge with opening angle ip = 20° for wall strength Ew = 0.9e. The 
filling transition at ksTj « 1.22£ is continuous. 


effects can lead to extra terms in the binding potential for wedge filling which may compete with those arising directly 
from the intermolecular forces. 

A very subtle question which we have not yet addressed concerns the precise nature of the change in order of the 
transition. In principle this may happen via one of two mechanisms: a tricritical point or a critical-end point. If 
there was a critical end-point then in the range of ip values where the filling transition is continuous there would still 
be a meta-stable low coverage state even at the filling temperature Tf. In Fig. 13 we show plots of the numerically 
determined grand potential, obtained via partial minimization, as a function of a constrained value of the thickness 
of liquid from the wedge apex, in the second-order filling regime. The left panel corresponds to T slightly below Tf 
while the right hand panel is slightly above Tf. Both graphs have a linear contribution oc (Tf — T)l^ proportional 
to the film thickness which changes sign at the filling temperature. This has a purely thermodynamic origin arising 
from the surface tensions and is responsible for the macroscopic prediction ([^. As the temperature is increased the 
location of the minimum smoothly increases and eventually disappears close to Tf when the linear term changes 
sign. It is clear there is no local minimum when T > Tf indicating that the change in order of the filling transition 
is via tricriticality. This means if we were to sit along the line of first-order filling transition temperature Tf and 
decrease the opening angle ip, the adsorption of the low coverage phase would diverge continuously as we approach 
the tricritical value of ip. However studying the nature of this divergence in more detail would be extremely difficult 
do to finite-size constraints. 

Finally, we mention that our density functional study is mean-field in nature and neglects long wavelength fluctuation 
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FIG. 11: Hysteresis in the adsorption as a function of temperatnre for a wedge, with opening angle ip = 120°, and wall strength 
= O.Se. 
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FIG. 12: (Color online) Coexisting equilibrium density profiles for a wedge with opening angle ip = 120° and wall strength 
Ew = O.Se at the filling temperature ksTf = 1.385£. 


effects associated with thermal wandering of the interface along the wedge. These certainly do not alter the location the 
filling boundary, 0[Tf) = (tt — ip)/2, which is determined by surface thermodynamics, nor the underlying mechanism 
for the change in the order of the phase transition, which depends on the competition between geometry and long- 
ranged intermolecular forces. The only influence of thermal fluctuations of any import concerns the roughness 
of the liquid-gas interface, which in the regime where the filling transition is second-order, is expected to diverge 
according to a universal power-law « (T/ This is not allowed for in mean-field density functional studies, 

which, as is well known, always yields an interfacial width of order the bulk correlation length. In the presence of 
long-ranged intermolecular forces however, the roughness ^_l, even allowing for interfacial wandering is always much 
less than the equilibrium film thickness of liquid adsorbed near the wedge apex, and mean-field predictions for all 
other quantities of interest should be reliable. 
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FIG. 13: Constrained grand potential as a function of the film thickness £u, for a right-angle wedge with the wall strength 
Em = 0.8e. For temperatures below the filling temperature ksTf = 1.378e (upper panel) the potential exhibits a single minimum 
which continuously shifts away from the wedge apex as the temperature increases. In the bottom panel the tail of the binding 
potential has a negative slope as required since T > Tf. There is no potential barrier between a local minimum and a global 
extremum, indicating the change in order is via tricritical point rather than a critical end-point. 


Acknowledgments 

A.M. acknowledges the support from the Czech Science Foundation, project 13-09914S. A.O.P. wishes to thank the 
support of the EPSRC UK for Grant No. EP/J009636/1 and EPSRC Mathematics Platform grant EP/I019111/1. 


^ R. Evans, Adv. Phys. 28, 143 (1979). 

^ J.S. Rowlinson and B. Widom, Molecular Theory of Capillarity, Clarendon, Oxford, (1982). 

® Fluid Interfacial Phenomena, edited by C.A. Croxton, John Wiley & Sons Ltd., New York (1986). 

^ S. Dietrich, in Phase Transitions and Critical Phenomena, edited by C. Domb and J. L. Lebowitz (Academic, New York, 
1988), Vol. 12. 

® Fundamentals of Inhomogeneous Fluids, edited by D. Henderson, Dekker, New York (1992). 

® R. Evans, J. Phys.: Condens. Matter 2, 8989 (1990). 

^ K. Binder, Annu. Rev. Mater. Res. 38, 123 (2008). 





13 


® D. Bonn, J. Eggers, J. Indekeu, J. Meunier, and E. Rolley, Rev. Mod. Phys. 81, 739 (2009). 

® W. F. Saam, J. Low Temp. Phys. 157, 77 (2009). 

D. E. Sullivan and M. M. Telo da Gama, in Fluid Interfacial Phenomena, edited by C. A. Croxton (Wiley, New York, 1985). 
M. Schick, in Liquids and Interfaces, edited by J. Chorvolin, J. F. Joanny, and J. Zinn-Justin (Elsevier, New York, 1990). 
R. Evans, U. Marini Bettolo Marconi, and P. Tarazona, J. Chem. Phys. 84, 2376 (1986). 

R. Evans and U. Marini Bettolo Marconi, J. Chem. Phys. 86, 7138 (1987). 

C. Rascon, A. O. Parry, N. B. Wilding, and R. Evans, Phys. Rev. Lett. 98, 226101 (2007). 

R. Roth and A. O. Parry, Mol. Phys. 109, 1159 (2011). 

A. Malijevsky and A. O. Parry, J. Phys.: Condens. Matter 26, 355003 (2014). 

A. O. Parry, A. Malijevsky, and C. Rascon, Phys. Rev. Lett. 113, 146101 (2014). 

P. Concus and R. Finn, Proc. Natl Acad. Sci. USA 63, 292 (1969). 

E. H. Hauge, Phys. Rev. A 46, 4994 (1992). 

K. Rejmer, S. Dietrich, and M. Napirkowski, Phys. Rev. E 60, 4027 (1999). 

A. O. Parry, C. Rascon, and A. J. Wood, Phys. Rev. Lett. 83, 5535 (1999). 

A. O. Parry, C. Rascon, and A. J. Wood, Phys. Rev, Lett. 85, 345 (2000). 

A. O. Parry, A. J. Wood, and C. Rascon, J. Phys.: Condens. Matter 13, 4591 (2001). 

D. B. Abraham and Maciolek, Phys. Rev. Lett. 89, 286101 (2002). 

D. B. Abraham, V. Mustonen, and A. J. Wood 63, 408 (2003). 

G. Delfino and A. Squarcini, Phys. Rev. Lett. 113, 066101 (2014). 

A. Milchev, M. Miiller, K. Binder, and D. P. Landau, Phys. Rev. Lett 90, 136101 (2003); Phys. Rev. E 68, 031601 (2003). 

M. Muller and K. Binder, J. Phys.: Condens. Matter 17, S333-61 (2005). 

N. R. Bernardino, A. O. Parry, and J. M. Romero-Enrique, J. Phys.: Condens. Matter 24, 182202 (2011). 

A. Malijevsky and A. O. Parry, Phys. Rev. Lett. 110, 166101 (2013). 

A. Malijevsky and A. O. Parry, J. Phys.: Condens. Matter 25, 305005 (2013). 

Y. Rosenfeld, Phys. Rev. Lett 63, 980 (1989). 



